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We study, by numerical simulations on a lattice, the behaviour of the gauge-invariant nonlocal quark conden- 
sates in the QCD vacuum both in the quenched approximation and with four flavours of dynamical staggered 
fermions. The correlation length of the condensate is determined to be roughly twice as big as in the case of the 
gluon field strength correlators. 
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1. INTRODUCTION 

The QCD sum rules approach to non- 
perturbative phenomena in strong interaction 
physics consists in evaluating the power correc- 
tions of the product of two hadronic currents via 
the Operator Product Expansion (OPE), putting 
the relevant, unknown information about the 
non-perturbative vacuum structure in the val- 
ues of local condensates ((g(O)g(O)), (G(O)G(O)}, 
etc.). However, it has been recognized in [|l] || 
that in many applications the x-distribution of 
the vacuum fields cannot be safely neglected: this 
leads to the study of gauge-invariant field corre- 
lators, or "nonlocal condensates" . 

The gauge-invariant two-point correlators of 
the gauge field strengths in the QCD vacuum, 
defined as: 



(x) = <0|Tr {G w (0)SG ua (x)S^} |0) , (1) 



where G MP = gT a G° p is the field-strength tensor 
and S = S(x,0) is the Schwinger phase operator, 
have been extensively studied on the lattice in 
the past M pf. The basic result is that the corre- 



lator T> fip ^ a {x), in the Euclidean theory, can be 
written as the sum of a perturbative-like term, 
behaving as l/|a;| 4 , and a non-perturbative part, 
which falls down exponentially 



^; P iW-exp(-| x|/A^), 

* Speaker at the conference. 



(2) 



with a correlation length A^ ~ 0.13 fm for the 
quenched SU{2) theory B, A^ ~ 0.22 fm for the 
quenched SU(3) theory [g§ and X A ~ 0.34 fm 
for full QCD (approaching the chiral limit) 0. 

Along the same line, in this paper we present a 
lattice determination of the quark-antiquark non- 
local condensates, defined as: 



Ci(x) 



/=i 



(3) 



where S(0, x) is the Schwinger line needed to 
make Ci{x) gauge-invariant and T l are the six- 
teen independent matrices of the Clifford's alge- 
bra acting on the Dirac indices a, b. The trace in 
(3) is taken with respect of the colour indices. 

The computations have been performed both in 
the quenched approximation and in full QCD us- 
ing four degenerate flavours of staggered fermions 
(whence the sum over the flavour index / in (3)) 
and the SU (3) Wilson action for the pure-gauge 
sector. We have adopted the same basic strategies 
and techniques already developed for the study of 
the gluonic correlators. 

In full QCD the nonlocal condensates have 
been measured on a 16 3 x 24 lattice at (3 — 
5.35 and two different values of the quark mass: 
a ■ m q = 0.01 and a ■ m q — 0.02. 

For the quenched case the measurement has 
been performed on a 16 4 lattice at j3 — 6.00, us- 



ing a valence quark mass a 



= 0.01, and at 
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Figure 1. The function a 3 Co(x) versus the dis- 
tance d, for the full-QCD case at (3 — 5.35 
and quark masses a ■ m q — 0.01 (circles) and 
a ■ m q — 0.02 (squares). The curves correspond 
to our best fits [Eq. (6)]. 
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Figure 2. The function a 3 C v (x) versus the dis- 
tance d, for the full-QCD case at f) = 5.35 
and quark masses a ■ m q — 0.01 (circles) and 
a-m q = 0.02 (squares). The curve corresponds to 
our best fit [Eq. (6)] to the data at a ■ m q — 0.01. 



(3 — 5.91 with a quark mass a ■ m q = 0.02. The j3 
values have been chosen in order to have the same 
physical scale as in full QCD at the corresponding 
quark masses, thus allowing a direct comparison 
between the quenched and the full theory. In this 
way it is possible to estimate the effects of the 
inclusion of dynamical fermions on the quanti- 
ties under study. In the quenched case two other 
measurements have been done at f3 = 6.00, using 
valence quark masses a ■ m q = 0.05 and 0.10. 

Further details about the computation and the 
results will soon be published |TlJ| . 

2. COMPUTATION AND RESULTS 

Making use of T,P invariance one can prove 
that all the correlators (3) vanish, except those 
with r* = 1 {"scalar nonlocal condensate") and 
with r* = 7^ and \x in the direction of x 
( "longitudinal-vector nonlocal condensate"): 

4 

Co(M) =-^(Tr[ ? t (0)5(0, x)q£(x)]) ; 
/=l 

4 

C M) = -o E^W^WW]} • (4) 

To discretize our operators on the lattice we 
have properly combined the staggered propaga- 
tors in order to build up the physical quark prop- 



agators. The Schwinger phase operator, i.e., the 
straight line of links connecting the point to the 
point x, has then been put in. In this way we con- 
struct two lattice operators C^d-a) and C^(d-a) 
(x — d ■ a, where d is the distance in units of lat- 
tice spacings a), which are proportional, in the 
naive continuum limit, to Co(d ■ a) and C v (d ■ a) 
respectively: 

C^d-a) - a 3 C (d-a) + O(a 4 ) , 

a— *Q 

CHd-a) - a 3 C v (d-a) + 0(a 4 ) . (5) 

Higher orders in a in (5) as well as possible mul- 
tiplicative renormalizations are removed by cool- 
ing the quantum fluctuations at the scale of the 
lattice spacing, as explained in Refs. g-|9j. This 
removal will show up as a plateau in the depen- 
dence of the correlators on the number of steps of 
the cooling procedure: our data are the values of 
the correlators at the plateaux. 

The nonlocal condensates have been measured 
at distances d — 2,4,6,8. In Figs. 1 and 2 we 
display the results for a 3 Co(d ■ a) and a 3 C v (d ■ a) 
respectively obtained in full QCD. 

We have tried a best fit to the data with the 
following functions: 

C {x) = A cxp(-^ x) + — ; 

x z 

C v (x) = A v x 3 exp(~fi v x) H j. (6) 
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Table 1 

Results for the correlation length Ao- Reported errors refer only to our determination and do not include 
the uncertainty on the physical scale ("/" stands for "full-QCD", while u q" stands for "quenched"). 



f3, theory a ■ m q Xq (fm) 



5.35, f 


0.01 


63 +U 2i 


5.35, f 


0.02 


46+ 09 


6.00, q 


0.01 


64+ ' 22 


5.91, q 


0.02 


46+ 006 


6.00, q 


0.05 


0.30(2) 


6.00, q 


0.10 


0.187(3) 



The form of the perturbative-like terms in Eq. 
(6) (i.e., Bo/x 2 for the scalar condensate and 
B v /x 3 for the vector condensate) is that obtained 
in the leading order in perturbation theory, in the 
chiral limit m q — * 0. Further details, as well as 
a remark about the reliability of the results ob- 
tained for the vector nonlocal condensate, will be 
soon published |ll]] . 

3. DISCUSSION 

The quantity of greatest physical interest which 
can be extracted from our lattice determinations 
is the correlation length Ao = 1 //j,q of the scalar 
quark-antiquark nonlocal condensate. It plays a 
relevant role in many applications of QCD sum 
rules, especially for studying the pion form factors 
and the pion wave functions In Table I we 

report the values obtained for Ao in all the cases 
examined. The physical scale in the full QCD 
case has been set by measuring and m p on 
our configurations, while for the quenched case 
we have used results reported in Ref. (!(]] . 

At a ■ m q — 0.01 the value is roughly twice as 
big as the value for the correlation length Xa of 
the gluon field strength at the same quark mass 
. As in the gluon case , the fermionic corre- 
lation length appears to decrease when increasing 
the quark mass. However, the full-QCD and the 
quenched values of Ao are nearly the same, when 
compared at the same quark mass, thus suggest- 
ing that this quantity is not much influenced by 
the introduction of dynamical quark loops. 

Finally, using the values of the pion mass m T 
measured on our configurations in the full-QCD 
case we find that m^Ao = 1.6(4) for a ■ m q = 0.01 



and TOttAo = 1-4(2) for a ■ m q — 0.02. There- 
fore, within the errors, the inverse of the scalar 
correlation length, fi = 1/Ao, turns out to be 
proportional to the pion mass m w . 
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